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Abstract 

The effect of quantum lattice fluctuations on the optical-absorption spectra in 
the ground state of halogen-bridged mixed-valence transition-metal linear-chain com- 
plexes is studied by using a one-dimensional extended Peierls-Hubbard model. The 
nonadiabatic effects due to finite phonon frequency uj n > are treated through an 
energy-dependent electron-phonon scattering function 8(k',k) introduced by means of 
an unitary transformation. The calculated optical-absorption spectra do not have the 
inverse-square-root singularity, but they have a peak above the gap edge and there ex- 
ists a significant tail below the peak, which are consistent with the optical-absorption 
coefficient or the optical conductivity measurements of this material. 
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1 Introduction 



Halogen-bridged mixed- valence transition-metal complexes (HMMC's or MX chain), as one 
of typical materials with a quasi-one-dimensional charge- density- wave (CDW) state, have 
been the subject of intense studies in recent years for both theoretists and experimentalists, 
because of the intrinsically interesting properties and the important technological applica- 
tions of the materials p], [3|]. Each chain of the materials is composed of transition-metal 
ions M +3 (=Pt +3 , Pd +3 , Ni +3 ) bridged by halogen ions X~ (=C1~, Br~, I - ), with ligands 
attached to the metals and in some cases counterions between the chains to maintain charge 
neutrality M +3 has an unpaired electron in its d z 2 orbital ("z is parallel to the chain), and 
this orbital makes an energy band through the supertransfer between neighbouring two d z i 
orbitals. This supertransfer comes from the hybridization between the d z 2 orbital and the 
p z orbital of X. Because of the electron-phonon (e-ph) coupling, the charge transfer occurs 
between neighbouring two M's so as to give the mixed- valence state. This is nothing but the 
CDW state with twice the period of the original lattice. The MX chain materials constitute 
a rapidly expending class of near single-crystal quasi-one-dimensional systems with phases 
of ground states ranging from strong CDW to spin- density- wave or spin-Peierls [Q . 

Many previous investigations of MX chain materials used the adiabatic approximation 
for the phonons. The adiabatic approximation means neglect of the dynamic energy of halo- 
gen ions, that is, neglect of the effect of lattice vibration. In nonadiabatic case, the phonon 
frequency u n is finite. Generally speaking, the nonadiabatic effect will suppress the state 
order parameters of the system || [|. As far as the optical-absorption coefficient is con- 
cerned, the results of adiabatic approximation methods have inverse-square-root singularity 
at the gap edge. By considering the nonadiabatic effect, the inverse-square-root singularity 
at the gap edge may disappear . The relationship between the change of optic-absorption 
shape and the phonon frequency in the range from u n = to u> w — > oo should be studied 
for understanding the physical properties of e-ph interactions of MX chain in nonadiabatic 
case. 
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The one- dimensional one-band model for MX chain complexes can be written as [§, ^, |Hj 

H = £ p f + ~ 51 *o(c!',* c »+i,« + c l+i,s c i,s) + a ( u i - u i+M,s c i,s + H c (1) 

where 



, v2M 1 2 



H C = UJ2 HtHl + V n i,s n i+i,s'- (2) 

Here c| a and q )S are the creation and annihilation operators of electrons sitting at site / 
with spin s, U\ (it's conjugated momentum is pi) is the displacement of the Zth ion, to is the 
supertransfer energy of electrons between neighbouring two orbitals, a is the e-ph coupling 
constant, K is the elastic constant and M the mass of ions. U and V correspond to the 
on-site and the nearest-neighbour Coulomb energies, respectively (throughout this paper, we 
set h = k B = 1). 

2 Theoretical analysis 

In the Hamiltonian ([!]), the operators of lattice modes, U[ and pi, can be expanded by the 
phonon creation and annihilation operators and b q , and after Fourier transformation, the 
H reads 

H = H° + H ep + H c , (3) 

where 

H ° = E w *( 6 & + ^+£e*4^« (4) 

q v z/ k,s 

H e P = Y,-m9{q){K + h - q ) c i+ q , s Ck,s- (5) 

q,k,s V ^ 

Here, = — 2to cos k is the bare band function, and iV is the total number of sites. The 



coupling function g(q) = ay l/(2Mov)[l — exp(ig)] and the dispersionless phonon frequency 

In order to take into account the fermion-phonon correlation, an unitary transformation 
is applied to H [[□], 12] , 

H' = exp(S)H exp(-S), (6) 



where the generator S is 



N q , k , s w w 



(7) 



Here we introduce a function £(£;', fc) which is a function of the energies of the incoming and 
outgoing fermions in the fermion-phonon scattering process and will be determined later. 
The unitary transformation can proceed order by order, 

H' = H° + H ep + H C + [S, H°] + [S, H ep ] + [S, H c ] + i [S, [S, H ]] + ± [S, [S, H c \] +0(a 3 ). (8) 

The first-order terms in H' are 

H ep +[S,H°] = lj29(q)(bl q + b q )cl +q ,sCk,s 

q,k,s 

7= E 9 (q)H k + 9, k )( b -g + b q) c l+ q ,s c k,s 



q,k,s 



1 ^ ^JL^fc + ?j _ bq ^ €k _ e k+q )c{ + c K 

q,k,s " 



(9) 



Note that the ground state \g ) of H°, the non-interacting system, is a direct product of a 
filled fermi-sea \FS) and a phonon vacuum state |p/i, ) fl3|| : 



\g ) = \FS)\ph,0) 



(10) 



Applying the first-order terms on \go ) we get 
+ tf°])|so } = E 9(q)blA +q ,sCk 

q ,k,s 



l-5(k + q,k)[l- 6k 6k+q ) 

/ 



l^o >, (11) 



since ) = 0. As the band is half-filled, the Fermi energy ep = 0. Thus c^g^IFS 1 ) ^ 

only if e k+q > and < 0. So, we have 



[H ep +[S,H°})\g ) = 0, 



(12) 



if we choose 

5(k + q,k) = 1/(1 + |e fc+9 - e k \/u v ). (13) 

This is nothing but making the matrix element of H ep + [S, H°] between \g ) and the 

lowest-lying excited states vanishing. Thus the first-order terms, which are not exactly 



canceled after the transformation, are related to the higher-lying excited states and should 



be irrelevant under renormalization |13fl . The second-order terms in H' can be collected as 
follows: 

[S,H ep ) + ~[S, [S,H ]} 
= ^EE ^^-6(k + q, k)[2 - S(k> + q>, k'))(bl ql + b q ,)(bl q - b q ) 

Zly q,k,sq',k' W * 

X ( c k+q,s C k',s8k,k'+q' ~ c k'+q' ,s C k,s$k' ,k+q) 

+4 E E 9 -^^S(k + q, k)5(k> + q', k'){e k , - e h , +q ,){bl q , - b q ,){bl q - b q ) 

ZJV q,k,sq>,k> 7r 
X { C k+q,s C k',s5k,k>+q> ~ C k'+q',s C k,s^k',k+q) 

~EE 9iqM ~ q) Kk + q, AO [2 - S(k' - q, k'U^A'-q^,,. (14) 

ZJV q,k,sk',s' ^ 

8k',k+q is the Kronecker 6 symbol. All terms of higher order than a 2 will be omitted in 
the following treatment. Note that the first order terms of transformed H c have not been 
included in the process of determining function S(k', k) (from Eq. (9) to Eq. (13)) and the 
second order terms of that have not been included in Eq. (14). The reason for doing so will 
be discussed later. 

Then we make a displacement transformation to H' to take into account the static 
phonon-staggered ordering 0, 

H = exp(R)H' exp(-R). (15) 

Here, 

R = -Y,{-^^^{b\-b l l (16) 

and exp(i?) is a displacement operator. 

If the ground state of if is \g), then the ground state of if is \g')\ \g) — exp(— S) exp(— R)\g'). 
We assume that for \g'), the fermions and phonons can be decoupled: \g') ~ \ fe) \ph, 0), where 
\fe) is the ground state for fermions. After averaging H over the phonon vacuum state, we 
get an effective Hamiltonian for the fermions, 

H cS = (ph,0\H\ph,0 



fc,s k>0,s 

4EE( + 9, fc) [2 - - g, fc')] - Vtay ) 

iV q,k,sk',s' I ^ J 

x 4+ g , s c fe, s 4'- g , S ' c fc', S '- ( 17 ) 



Where 



Vo = i ^ g ^ (fc + g,fc)[2-^-g,fcO]- (20) 



= efc-^E ^' ^ (18) 
A (k) = 2au [l-S(k-Tr,k)}, (19) 

k,k',q 

Now, the total Hamiltonian can be divided as H = H + Hi , where Hi includes the terms 
that are zero after being averaged over the phonon vacuum state, and 

H = J2^(b% + ^)+H cS . (21) 

After the unitary transformations, the Coulomb interaction H c becomes if c +higher order 
terms. The contributions of these higher order terms to H Q are zero under the mean-field 
approximation. So, in the process of determining function 5(k',k)), which only concerns 
the first order terms of transformed Hamiltonian, and also in Eq. (14) and Eq. (17), we 
have not written out the higher order terms of transformed H c . By using the mean-field 
approximation to the terms of four-fermion interaction in H e n, we have 

H = \KNul (b\b q + 1) + E^4,s c M + E A M- n , s Ck,s + 4,-Cfc-*,.). (22) 

A q \ £ ' k,s k>0,s 

In the theoretical analysis, we have distinguished different physical processes and taken into 
account the fact that only the Umklapp scattering terms affect the gap, and the forward and 
backward scattering terms contribute nothing to the " charge" gap [|14], |T^j . The renormalized 
band function and the gap function are 

Ek = E (k)-^Y,{^^ 2 (^n k ^m~S(k\k)}~Vcos(k'-k)\^-, (23) 

k 1 \ ) " 

A k = 2au [c--8(k-7r,k)}. (24) 



Where 

A, 



1 + ^ 

iV fc>0 



2V-\u-V + ^6(k-n,k) 

Z K 



2auoWk ' 



(25) 



and Wk — y -Ef + A| is the fermionic spectrum of elementary excitations in the dimerized 
state. 

We find by means of the variational principle that the equation to determine dimerized 
lattice displacement ordering parameter uq is 

^J^-^^^k (26) 



If uj w = we have 5(k', k) = 0. Eq. (|26j ) becomes the same as that in the adiabatic theory. 
In our theory 5(k — n, k) has a sharp peak at the Fermi point and, since 

u;^ + 4to| cos k\ 



the logarithmic singularity in the integration of Eq. ( [ZED i n the adiabatic case is removed 
by the factor 1 — S(k — n,k) as long as the ratio u^/to is finite. For the gapped state 
(insulating state) we believe that the mean-field approximation is still an effective one. 
The reason is as follows: The quadratic part of H e g (which can be exactly diagonalized) 
contains a gap term (Ao(k) ^ 0) that is related to the Peierls dimerization and suppresses 
the quantum fluctuations. Furthermore, the logarithmic singularity (related to the infrared 
divergence of low dimensional system) has been removed by the factor 1 — 8{k — tt, k) = 
4t | cos k\/(u n + 4t | cosfc|). Comparing Eq. (|24] ) with that in the adiabatic case, A = 2cra , 
we have the gap in the nonadiabatic case, A = A(7r/2) = 2auo[c — 1/2], which is nonzero 
but smaller than the adiabatic gap. This is the true gap in the fermionic spectrum. 

3 Optical absorption 

The optical-absorption coefficient a{u) can be expressed by the retarded Green's function 
as follows: 

a(u) = — -lmK R (u), (28) 

71UJ 
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where K R is defined as 



K R (u) = -i f e-^dt(g\[j(0)j(t) - j(t)j(0))\g). (29) 

<J — oo 

Here, j is the current operator 0, 

j = -iet J2(4, s c W,s ~ 4+i,s c i,s) 

1,8 

= Jj2 sinkc k,s c k,s, (30) 

where J = 2eto, and j(t) = exp(iHt)j exp(-iHt) is the form of j in the Heisenberg rep- 
resentation. As R commutes with j, the unitary transformation of the current operator 
is 

e S je S = J + [S ^ ]+ 1 _ [Si[Sij]]+0{a 3 ) 



k>0,s 

+-7T7 Yl —( b -g ~ h i)K k + 1, k)[sink - sin(A; + 9)1(4+ c M - c l+ q ~n,s c k-n,s) 

J ^ g(q)g(q') (fct _ g _ b9)(p i _ bqi)§{k + g> k)[sink _ sin{k + g)] 

2iV k>0,q,q',s U * 

x[5(k + q + q',k + q){cl +q+ql>s c kjS - c\ +q+ql _^ s c k ^ s ) 

-S{k, k - q'){c\ +qs c k _ q ,, s - 4+ <z _ 7r , s Cfc- g '- 7r)S )]. (31) 

Because the averaging of Hi over the phonon vacuum state is zero, in the ground state at 
zero temperature, Hi can be neglected. By using the approximately decoupling \g') « | <7q) , 
the ground state of H , and H « H |13| in the calculation 



(g\mm\9) = ^|[ e (5+/?) je- (5+fl) ]e^[e^)je-^)]e-^|^) 

« W \[^je- 8 ]e^Vje- S ]e^y ), (32) 



we can get 

k» = J 2 E 



fc>0,s 



f 2a 2 1 A 2 

x ^ sin 2 k + — — — Et 1 - cos ( fc ' - ^^(^ fc ) sin M sin - sin k ) \ 7772 
{ KNu) n k , J W k 



a 2 J 2 _ 

+ KN 2 cu„ fc>0 ^ >0jS \u-u n -W k -W k , + m u + u n + W k + W k < - zO+ 
x [[1 - cos(At' - k)]8 2 (k',k)(sink' - sink) 2 (a k p k > + (3 k a k i) 2 
+ [1 + cos (A/ - k)]5 2 (k' - 7T, A;) (sin A;' + sin k) 2 (a k a k , + /3 k (3 k >) 2 } , (33) 



where a k = y (1 + E k /W k )/2 } and = y (1 — E k /W k )/2. Thus, we have the optical- 
absorption coefficient 

a{u) = 2j2{2S+1) J2^-™ k ) 

WiV fc>0 

f 2a 2 1 A 2 

x <^ sin 2 k H Vfl - cos(A; / - k)]5 2 (k' k) sin A;(sin kl - sin fc) \ — | 



+ £ 



2a 2 J 2 (2s + l) 

k>0,k>>0 

x |[1 - cos(A; / - k)}5 2 (k', k)(sink' - sink) 2 (a k l3 k > + (3 k a k f) 2 

+ [1 + cos(A/ - k)]5 2 (k' - tt, A;) (sin A;' + sin k) 2 (a k a k , + (3 k (3 k ,) 2 ) . (34) 

The relationships between the change of the optical-absorption shape and the different 
phonon frequencies are shown in Figure 1. The parameter values used are: a 2 / K = 0.5*o, 
U = l.Oto and V = 0.15* with to n = 0.004* , 0.008* , and 0.012* . One can see that the 
spectrum broadens but the peak height decreases as u n increases. The inverse-square-root 
singularity at the gap edge in the adiabatic case |T(| disappears because of the nonadiabatic 
effect. For comparison, the adiabatic (uv = 0) result is also shown which has an inverse- 
square-root singularity at gap edge. For finite u n , the singularity is absent and there is a 
significant tail below the peak. We note that in our theory, in mathematical viewpoint, the 
difference between the u n = and u n > cases mainly comes from the functional form of 
the gap [see Eq. (E3j)]. Comparing it with that in the adiabatic limit, one can see that the 
subgap states come from the quantum lattice fluctuations, i.e., the second term in the square 
bracket of Eq. ([§). 

Figure 2 shows the calculated optical-absorption coefficients a; (a;) versus the normalized 
photon frequency relations for u n = 0.008* , U = 1.0* and V = 0.15t m the cases of 
a 2 1 K = 0.3t , 0.4* , and 0.5* . One can see that as the e-ph coupling constant a 2 / K 
increases, the peak of optical-absorption spectrum moves to higher photon energy and the 

9 



dimerization gap becomes wider. 

Figure 3 (a) and Figure 3 (b) show the calculated optical-absorption coefficients a(ui) 
versus the normalized photon frequency relations for u n = 0.008to and a 2 / K = 0.5£o in the 
cases of different U and V values, respectively. One can see that the spectra are generally 
the same, except that the peak height decreases and shifts to higher photon energy when 
U decreases or V increases. The on-site and the nearest-neighbour Coulomb repulsions 
are to weaken and strengthen the lattice dimerization, respectively. In the MX chain, a 
charge-transfer (CT) excitation originates in the electron excitation over the CDW gap. If 
no definite fine structure can be observed in the absorption band, it cannot be clarified, only 
from the line shape of optical-absorption, whether this band is due to the CT exciton or the 
interband transition. Wada et al. 0] claimed that because of the electric field dependence 
of the modulation signal observed in the measurement of the electroreflectance spectra of 
the MX chain, it can be excluded that the CT excitation absorption band is an interband 
absorption. However, K. Iwano et al. fl8|] indicated that in the case of weak electron-hole 
binding, i.e., small V, the spectrum becomes the mixture of the exciton absorption and the 
interband ones. If the value of V is reduced toward zero, the exciton contribution then 
becomes more and more small and, finally, the exciton disappears and the spectrum is now 
given completely by the interband transitions. Our calculated result of the V = case is also 
shown in Fig. 3 (b) and this is the optical-absorption spectrum of the interband transitions. 

In Figure 4, we compare our calculation (solid line) with the observed optical absorption 



(solid circles) of MX chain complexes ||19| . Here we use the same input parameters (without 
adjusting them) U = l.Oto and V = 0.3to, as those of K. Iwano and K. Nasu 18|| . Other input 
parameters used are u n = O.Olio and a 2 / K = 0.2to- One can see that the agreement between 
the experiment result and our calculation is quite good. In experiments, the measurement 
of optical-absorption spectrum could be affected by various factors, such as the impurity in 
samples, the finite measurement temperature, and the finite measurement resolution (about 
±5% accuracy of the reflectivity measurement [fjl|]). These made the measured absorption 
spectrum to be broadened and led to the slight discrepancy between the experiment result 
and our calculation. 
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In Figure 5, we plot the calculated optical-absorption coefficients versus the rescaled 
photon frequency relations with U = l.Oto and V = 0.15£o for different u n and a 2 / K. Quite 
remarkably, when we scale the optical-absorption coefficient by the peak value, and scale 
the photon frequency by T, the half- width for the low photon frequency side of the peak, 
we find that the scaled curves for various u> n and a 2 /K values have an universal form below 
the peak photon frequency. In their studies of temperature dependence of subgap optical 
conductivity for KCP(Br) and trans-polyacetylene, J.W. Wilkins et al. [f2(]] discovered that 
the strong subgap tail of the conductivity had an universal scaling form. We use the same 
scaling method as that of J.W. Wilkins et a/., but to study the u n and a 2 / K dependence of 
optical-absorption coefficient. Our universal form is similar to that of J.W. Wilkins et al, 
which implies that the temperature effect and nonadiabatic effect on quasi-one- dimensional 
Peierls systems have similar characteristic properties. 

4 Conclusions 

The effects of quantum lattice fluctuations on the optical-absorption coefficient in the ground 
state of MX chain are studied by introducing an energy- dependent e-ph scattering function 
8(k',k) via an unitary transformation. The functional dependence of 5(k',k) on k! and k 
is determined by the vanishing of the second-order perturbation of first-order terms in the 
transformed Hamiltonian. By using the retarded Green's-function method we have shown 
that our theory gives a good description of the effects of quantum lattice fluctuations and 
change of various parameters on the optical-absorption. The inverse-square-root singularity 
of optical-absorption spectrum at the gap edge in the adiabatic case disappears because 
of the nonadiabatic effect. The calculated optical-absorption coefficient is consistent well 
with the experiment observation of optical conductivity measurements. The curve of subgap 
optical-absorption coefficient has an universal scaling form. 
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Figure Caption 



Fig. 1. The relationships between the change of optical- absorption shape and the dif- 
ferent phonon frequencies in the cases of a 2 / K = 0.5* , U = 1.0* and V = 0.15* with 
cov = 0.004* , 0. 008*o, and 0.012* . For comparison, the adiabatic (ui n = 0) result is also 
shown which has an inverse-square-root singularity at gap edge. 

Fig. 2. The calculated optical-absorption coefficients a(uj) versus the normalized photon 
frequency relations for u n = 0.008*o, U = 1.0*o and V = 0.15*o in the cases of a 2 / K = 0.3*o, 
0.4*0, and 0.5* . 

Fig. 3 (a). The calculated optical-absorption coefficients at(u) versus the normalized 
photon frequency relations for u n = 0.008to, u 2 / K = 0.5to and V = 0.15to i n the cases of 
U = 0.6t ,0.8t , and 1.0* . 

Fig. 3 (b). The calculated optical-absorption coefficients a(u) versus the normalized 
photon frequency relations for uo^ = 0.008*o, a 2 / K = 0.5*o and U = 1.0*o in the cases of 
V = 0, 0.15*o, 0.25*o, and 0.35* . 

Fig. 4. The rescaled optical-absorption coefficient a (a?) / 'a pea k versus the normalized 
photon energy u/upe^ relations, where a pea k and Wpeak are the optical-absorption coefficient 
and the photon energy at the peak, respectively. The solid circles are the result of experiment 
observation for [Pt(en) 2 ][PtX 2 (en) 2 ] (0104)4 (Ref. [fL9|1 , Fig.6). The solid line is our result by 
using ccv = 0.01*o, U = 1.0* , V = 0.3* , and a 2 / K = 0.2* . 

Fig. 5. Scaling plot of the calculated optical-absorption coefficients a(dj)/a peak , where 
the photon energy is scaled by T, the half-width for the low photon energy side of the peak. 
Note we claim only that the subgap(u; < Wpeak) optical-absorption coefficients scale. 
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